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COORDINATE AXES, LOCATION OF ORIGIN, AND REDUNDANCY 
FOR THE ONE AND TWO DIMENSIONAL DISCRETE FOURIER TRANSFORM 

FOR COMPLEX AND REAL DATA 

George E. I6up arid Juliette W. loup 
Department of Physic* and 
Geophysical Research Laboratory 
University of New Orleans, New Orleans, LA 701 AS 

For continuous Fourier transforms, the shift theorem will 
locate the origin anywhere. If 

*Cf<x>> - F <s> , 

then 

^Cf(x-a>> ■ exp (-iSiras) F<s) 

for ’forward transforms. For inverse transforms, 

?~1<F <s-a) )■ * exp< + iSTras) f(x) 

The discrete Fourier transform (DFT) has special 
considerations: (1) there are only a finite number of sample 

points, and <E> there is replication in both domains. 

COMPLEX ONE-DIMENSIONAL DATA . 

Consider a function domain array ax, where x gives the 
coordinate location of each of the complex elements of ax* The 
transform of ax is given by the complex array Aj. Here j 
specifies the frequency number defined by jAs = s. For Ax = 1, 
As is numerically equal to 1/N, with N the number of points in 
a x . Its units are the reciprocal units of x. In general 
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As a l/<NAx> 

Data arrays often have no particular origin associated with 
their elements. Those that do are commonly one of two forms: 
Form 0a: 

a 0 a l as • • • a N-£ a N-l 
which h 'is N points, and 
Form las 

a-M a-M+1 ••• a -£ a-l a 0 a i as • • • a M-l a M 

which has <£M+i) points. The standard convention for the DFT 

is 

N-l 

Aj ® £ ak exp (-i£trkj/N) , 

k=0 

where k gives the position of a*< as measured from the left of 
the array (with the first element having k ** 0) as stored by 
existing fast Fourier transform (FFT) programs. The form of 
the exponent is such that k = 0 is assumed to be the label of 
the origin element. Data in Form 0a can be directly 

transformed, but for data in Form la, something must be done to 
have the transform correctly represent the data with their 
given origin. As an alternative to having different FFT 
algorithms for different data origins, the data may be modified 
to use existing codes. 

The required modification of the data is a rearrangement 
as follows: 

Form 3a x 

*0 *1 *2 . . . a M - 1 atv] a-|r |+ 1 ... a-£ a - 1 

for N ■ £M+1, or 
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Form 4a: 

30 3 J i • ■ 3|y(— 1 *"M 3— (vj4* J » • • 3-— g 1 

\ 

for N “ £M, with am» ~M * a|v| + a-M- N is the number of points 
to be transformed. 

Because the DFT corresponds to replicated data, the above 
modifications describe a shift of the data to the left by M 
elements, or of the origin to the right by M elements for N 
even, or by M+l for N odd. If N is even, another form is 
needed to correspond to Form la. The first element of Form la 
is replaced by aty , -m and the last element dropped to give Form 
£a. 

Form £a: 

am, -M a-M+i ... a-£ a-i a 0 ai ag ... aft-i 

Another approach consists of taking the DFT of the 
unmodified data directly. The DFT definition is then used in 
it* standard form and the DFT is corrected, rather than 
rearranging the data before the DFT is calculated. The DFT of 
the unshifted data can be changed into that of the shifted data 
by a phase multiplication determined by the Shift Theorem 
(Bracewell, 1978). Since the origin is assumed to be at the 
left-most element for the standard DFT algorithms, the a* array 
must be shifted to the left so that a 0 will fall on the origin. 
For the forward DFT (-i transform), the phase multiplier 
corresponding to this shift will be exp ( + i£irsMAx) . The 
multiplication of the DFT of data in Form la or £a by this 
factor will result in a DFT corresponding to an origin 
correctly located at ag). 

The transform of sampled daita is replicated, which means 
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that the DFT is periodically repeated. The period is 1/Ax, 
which is often assumed to be 1. Since a periodic function is 
completely specified by one period (called an "island" by 
Bracewell, or principal interval, or base band), the discrete 
transform is given as one period of the periodic result. The 
base island of the transform, centered on the origin, is the 
part of the transform generally used for conceptual 
manipulations. Working with this island has the advantage that 
any frequency-dependent filtering will be simplified since the 
s values (frequency values) are simply associated with the 
elements of the transform arrays 
Form 1R* 

R— M ... R— •£ R— 1 Ri0 Ri R£ ... R|*i— 1 

s= -MAs (-M+1 ) As. . . -2As -As 0 As iSAs. . . (M-l) As MAs 

for N odd. Unfortunately, the convention for the DFT is not 

the above arrangement, but is 

Form 3fls 

R0 Ri Rg . . . Rm~ l R|v| R~m+ 1 • . « R-£ R- 1 

s» 0 As £As. .. (M-l ) As MAs (-M+1 > As. . . -2As -As 

for N » 2M+1, or 

Form 4Rt 

R0 Rl R£ . . . Rm-1 Rm» -M R-M+l fl-2 0-1 

s* 0 As 2As. . . (M-l) As -MAs (-M+1 ) As. . . -2As -As 

for N => 2M, with flM* -M 5 RM + A-M* This arrangement may be 

thought of as consisting of the origin and right-hand half of 

the base island and the left-hand half of the first replication 

to the right. When filtering the transform in this form, care 

must be taken to associate the correct s value with each 

element. Rs an alternative to working with this DFT result, 
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the elements of the DFT may be rearranged into Form 1A for N 
odd or into 
Form 2A: 

Am, — M A— (*i-f 1 • ■ • A— 2 A— 1 A0 Ai Ag • ■ * Am—1 

for N even. This may be accomplished by rearranging the DFT 
elements after they are obtained, or by phase multiplying the 
data before the transform is taken. Again the shift theorem is 
used to determine the phase multiplication. The phase 
multiplication for each a* will be exp <-i£ 7 rM (As) x > if the -i 
transform is used for the DFT. 

In the transform domain, the phase multiplication needed 
to correct for Form 2a with N even <M = N/2> is 

eiSirWsAx k *i2ir<N/2) JAsAx s e i£ir(hi/2) JAx/ (NAx) 

This becomes 

eiirJ = (-l)J 

Care must be taken that j measure the displacement in terms of 
positive or negative position number of the elements from the §, 
origin of the transform. For <-l)J, positive or negative makes 
no difference, but for other cases it will be important. The 
phase multiplier for N odd also simplifies, although not as 
much. This phase multiplication corrects fo.r a data origin 
given by Form 2a. 

For the function domain, the phase multiplication needed 
to have the transform be in Form 2A, N even, is 
«-i£7r<MAs)kAx * e -i2ir(N/£) <l/NAx)kAx = e -i-rrk = (-i)k 

Again k measures the positive or negative displacement relat ive 
to the x origin . The simplification for N odd is similar but 
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not as great. This phase multiplication causes the transform 
to have Form £0 rather than Form 4A. 

For Form 4A or SA, for filtering purposes the element 
Am, is associated with s »* +MAs and s * -MAs since ft Mi ~M * ^M 
+ A~m- If the filter is H(s>, then the filtered result at this 
location should be 

H (MAs) Am. -m + H ( -MAs ) Am. ~M 
£ £ 

if the a>< elements are real, since the transform is Hermit i -an 
(conjugate even or symmetric, Bracewell, 1978). For complex 
ax, no rule can be given. Setting the filtered result to zero 
or to the result which holds for real data are two possible 
approaches. 


TWO-DIMENSIONAL DATA 

Consider Form 0a. Notice that the as, 0 element is in the 

* 

upper left corner, and subscripts increase down and to the 
right. This arrangement is chosen for two reasons? first, this 
is the usual way matrix elements are subscripted, and second, 
this is the usual order in which image data are stored in the 
computer. The subscripts can be used to denote the x and y 
axes of the two-dimensional transform. There, are two choices: 
ax,y : the x axis is vertically down and the y axis horizontal 

to the right, with the z axis out of the paper? or ay, x * the x 
axis horizontal to the right and the y axis vertically down, 
with the z-axis into the paper. The latter convention is 
chosen here because the x and y axes are generally taken to be 
horizontal and vertical, respectively, - in two-dimensional 
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coordinate systems. Note that the limits on >< and y are 
0 & x i Ng-1 or -Mg & x & Mg and 0 i y £ Ni~l or -Mi 1 y i Mi. 
Forms 0a or 3a are assumed by standard FFT algorithms. In the 
notation of Bracewell the exponential in the two-dimensional 
Fourier transform is given by exp (-Srri <xu+yv) ) . If x is chosen 
to be horizontal? then the u axis is horizontal j and if y is 
chosen to be vertical, then the v axis is vertical. 


Form 0a s 




T a0,0 
1 

*0, l 

a0i 2 • • • 

a0, NS-1 1 

1 

1 «lt0 
1 

aii 1 

a 1 , g ... 

«1,N8-1 * 

1 

1 ag,0 
1 • 

ag, 1 

ag, £ 

ag, NS-1 f 

a 1 

1 * 

I • 

1 *N1~2*0 
1 

a|\|l~£i 1 

a Nl-Si £ • ■ • 

• 1 
a Nl-£i NS-1 1 

1 

L »Nl-li 0 

*Ni-l» 1 

a|Ml-li£ • • ■ 

»Nl-li NS-1 J 

Form lai 

Mi, Mg 

odd 



r a-Ml 4 -MS 

a-Mli -MS+1 

a-Ml, -1 

a-Mli 0 
1 

a-Mli 1 • ■ • 

a-Ml 1 MS-1 

ja-Ml + la -MS 
1 * 

' • 

1 a 

a— M 1 + 1 , "-MS+ 1 • ■ 

III 

• a-Mi+li -1 

a-Ml+li 0 
1 
1 
1 

a-Ml + li 1 * « • 
IV 

a-Ml+li MS-1 

1 *-l » -MS 

a— 1, —MS +1 ■ ■ • 

a-1,-1 

1 

a-1,0 

1 

a-l, 1 ■ - • 

a-l, MS-1 

) -a0, —MS 

-a0, -MS +1 

a©, -1 

"*0, 0— — 
| 

-a0, 1 

a©, MS -1 

1 ai, -MS 

i ■ 

V • 

1 a 

ai, -MS +1 ••• 
II 

• 1, -1 

*1,0 

1 

1 

1 

a l , 1* • ■ • 

I 

«li MS -1 

I ami-1, -MS 

aMl-li -MS +1 ■ ■ 

■ »M 1 - 1 1 - 1 

1 

aMl-li 0 
| 

*Ml-li 1 ■ • ■ 

«M 1 — li MS -1 

L «M 1 t -MS 

*M 1 1 -MS +1 • • ■ 

»Mli -1 

*Mli 0 

*Mli 1 • • • 

aMliMS-1 


»-MliM£ 1 

a-Ml+li MS I 

• I 

a 1 

« 1 

a- 1» MS I 

*0» MS I 

*1 1 MS I 

■ I 

• 1 

e I 

*Ml-li MS I 

*Mli MS J 
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Form 2a s Mj, Ms «v®n 


*±M1§ ±M2 

*±M1 i -M2-H • • • 

a±Ml i -1 

1 *±M1 1 0 
1 

a ±Mii l « • ■ 

*±Mli M2-1 

1 

I 

*~Ml+li ±M2 

m 

a 

* 

*~Ml+li -M2+1 
III 

a-M 1 + 1 i "* 1 

1 *-Ml + li0 

1 

1 

*~Ml+li 1 • ■ • 
IV 

*-Ml+l i M2-1 

a 

a 

a 

1 

1 

§ J : 

*-li ±M2 

*— li -M2+ 1 • • « 


1 

1 a ~l»0 

1 

a— i » i » • » 

a -li M2-1 

! ( 
i t 
i 


— _ — 

- — i - 

| “ 

— 

— — — 

*0i ±M2 

*01 -M£+ 1 • • * 

*0i-i 

1 a®, ® 
1 

*01 1 • * * 

*0i M2-1 

! ! 

*li ±M2 

a 

a 

a 

*1 i -M2+1 • • • 

II 

*li-l 

1 a 1 1 0 
1 
1 
1 

a l* 1 • • • 
I 

*1 i M2-1 

a 

■a 

a 

i ■ 

i i; 

1 l 

1 j- 

*Ml~li ±M2 

a Ml“i i “M2’fl ••• 

*Ml-li -1 

I 

• awi-1,0 

a Ml-li 1 • ■ • 

*Ml~li M2-1 

J ! : 


Form 3a t Mj, Ms odd 


r * 0.0 

a 0i 1 * * * 

*01 M2-1 

*0i M2 

*0i -M2 

*0i — M2+1 " • • 

*01-1 

1 

• a 1 1 0 
1 a 

1 * * 

| a 

*1 i 1 • ■ ■ 

a l i M2-1 
I 

*li M2 

*li -M2 

* 1 i — M2+ 1 • • • 
II 

*li-l 

a 

a 

a 


1 a Ml — 1 1 0 

a 

3 

h* 

1 

« 

K 

» 

a 

■ *Ml-li M2-1 

*Ml-li M2 

*Ml-li-M2 

*M1— li -M2+1* • * 

*Ml-li -1 


1 *Mli0 

a Ml i 1 • • • 

*M1 * M2-1 

aMl i M2 

*Mli -M2 

*Mi » -M2+1 • » • 

*Mli-l 


1 *-Mli0 

a — Ml i 1 * * • 

*-Ml iMS-1 

*-MliM2 

*-Ml i -M2 

*— Mli — M2+1 • ■ • 

*-M 1 1 — 1 


I*-M1 + 1i 0 
i ■ 

1 ■ 

1 a 

a-Ml + li 1 • ■ 

• *-Ml+l i M2-1 
IV 

*-Ml+liM2 

*-Ml+li-M2 

a-Ml+li -M2+1* ■ 
III 

• *— Ml + 1 1 — 

a 

a 

a 

1 1 

L *-li 0 

*-li 1 ■ ■ • 

*-li M2-1 

*-l i M2 

*-li -M2 

*— 1 i — M2+1 * • • 

*-!• -1 

J 


a 


Form 4a j Ml* Mg even 


r *0, 0 a®, 1 . . . a«,M2»l • 

I 1 

I *1*0 *1*1 **• *1*M£”1 I 

I I I 

I . I I 

| a ) 

I I 

I *M1-1 * 0 a Ml-l* 1 ••* a Ml-l* M£-l I 

I I 

I a±Ml*0 * ±M 1 * 1 ••• a±Ml* M£-l I 

I I 

I a-Ml + 1,8 a-Mi + i,i.., a-Ml+liME-l I 

! a I 

I . IV i 

I • * 

l l 

La- 1,0 a-i* i ... a-l,M£-l 1 

with «|<f ±M£ - *k*M£ + «k»-M£ 
and a±Ml»±M£ “ »M1*M£ + *-Ml»M£ 


*0* ±M£ ®0*-M2+l *3»-l 1 

* 1 * ±M£ *1 * -M£+l *1.-1 < 

• I 

II . I 

i 1 

aMl-1* ±M£ «M1-1* -MS-fl «Ml-l*-i I 

a±Ml* ±M£ **Ml*~MS+i ••• *±M1* -1 I 

a~Ml+l* ±M£ *-Mi+l* ~ME+1 a— mi+i,~i I 

III . ! 

• I 

a-1* ±M£ «-l»-MS+l ••• a-l*-l 

a±Ml* I - *M1» I + *-Ml* I 
*M 1 * 1 -M£ + *-Ml* -ME ■ 


Tha quadrant numbers are shown as Roman numerals for 
Forms la, £a, 3a, and 4a. The quadrant numbering system 
depends upon the choice of the x-y coordinate system discussed 
above. Elements a©, i and ai ,0 are o& the axes. They are 
between quadrants and are not actually in. them. The rearranged 
sections from Form la to Form 3a are not all of the same size; 
they do not have the same numbers of columns and rows. Note 
that for N even, Form £a, the origin is not at the geometric 
center of the matrix, bu«i is just to the right of and just 
below the center. 

One may rearrange Form la into Form 3a or Form £a into 
Form 4a and then take the standard OFT. The result will 
represent the .data with the correct origin of Form la or Form 
£a. Alternat ively, one may take the transform of Form la as it 


A < 
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is for Nj, Njg odd (or Form S* for Ni» Ns eve n), and then phase 
multiply the transform to have it correctly represent the 
origin as in Form la or Form £a. 


In one dimension, the data are rwplicated and aliasing 
may occur. Two-dimensional data are also replicated and 
similar aliasing problems may be present. 
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II 
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in 

IV 

III 





Form la > 

VB 

1 

1 I 

li 

II i 

. r 


3 a ^ .. -- 




hi 

1 IV 

1 

HI | 





• 








Form la has M o" marking its center, and "x" marks the center of 
Form 3a. 


The phase multiplication in the transform domain to get 
the transform of data with the origin in the center will be 
exp(+£iriuMsAx) exp <+£TrivMiAy> , 

where the uv origin is in the center or the upper left corner 
of the matrix, depending on whether the transform has been 
rearranged or not. Since for Ni and Ns even, Ms * Ns/£, Mi * 

Nj/£, and since Au ■ l/(NsAx>, and Av * l/(NiAy), then the 
exponentials above become 

exp(+£iri <NsAx/£> (s/NsAx) > exp(+£iri (NiAy/£) (r/NlAy) ) 

■ exp(iirs) exp(iirr) - <-!)* <-i>r » <-l>r+s , 
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where r and a should be measured from the origin. 

The normal way to think of the transform is Form 1ft, 
which may be obtained from ths diagram for Form la by 
substituting capital latter ft for small a. However, the 
standard DFT calculated is Form 3ft for N odd and Form 4ft for N 
even. 

For Form 4ft, the dc level, ft«, 0, is in the upper left 
corner, and lower frequencies, i.e., those of small lul,!vl 
values, are found in the corners of the matrix. The highest 



Therefore, for two-dimensional low-pass filters, the center 


region must be blocked out. Note that the filter covers a two- 
dimensional area. 

If the transform origin is to be in the "center", 
rearrange the transform after taking it, or phase multiply the 
data before taking the transforms 
e x p < -2iri M£ AumAx ) ex p < -Siri M i Avn Ay ) 

* exp (-2-iriNsmAx/ <£NsAx) ) exp (-SiriNinAy/ <£NiAy> ) 
m <-i)m+n for Ni, Ns even 

The transform will be in Form 1ft for Nj, f Ns odd, or Form £ft for 
Ni, Ns even, for either method. 

Mith currently defined FFT computer algorithms, the 
initial data are assumed to be in Form 0a or Form 3a for Ni, Ns 
odd and Form 0a or Form 4a for Ni, Ns even. The distinction 
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between Forms Ba and 3m or 4a does not matter if the data are 
never considered mm giving the "frequency content" of the 
f »onward transform. 

C1NE-D I MENS I ONAL REAL DATA 

Whether calculating the transform of two real arrays with 
one complex transform, or the transform of on* N-point real 
array with on* (N/S) -point complex transform (Brigham, 1974) , 
or the ordinary complex transform for real data, computer 
storage requirements can be rediiced by making use of the 
symmetries of the transform for real data. The function domain 
forms do not change. The transform domain forms can be 
modified as followsi 
Form lAi 

A— M O-M+i * * * A-2 A— i A© 

♦ 

for N odd, with <<N~1>/S> + 1 complex term*. It seems that 
there are £C<(N-l)/2) + 13 ■ N + 1 independent pieces of data, 
which is one too many; but A 0 is real, so that there are 
actually only N independent values as there should be. To fill 
in the rest of the terms, make use of the Hermit ian propertyi 

# 

At ■ A~i 

Form 3A» 

A© Al ■ AM-1 am 

for N odd. There are N numbers represented as in Form 1A. 
Again use the Hermit ian property to fill in the missing termsi 

# 

A-i - Ai 
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Form 4ft: 


fi® fti fi£ • • • fiM— i Am* — m 

for N even, where both 0® and ftM, ~M are real. (This again 
insures that there are £C(N/£) 13 -■ £ ® N independent pieces 

of information.) 

Form £ft: 

fiM* — M fi-M+1 ... fi-£ fi-1 0® 

For real data one may choose to work only with positive 
frequencies, Form 3ft or Form 4ft. 

CfiLCLJLftTION OF THE TWO-D I MENS I QNftL DFT 

Recall that the definition of the two-dimensional DFT is 
Ni-1 N£-l 

firs “ £ £ am exp-C-iEsT (rk/Ni ) ■+• <s I/Ng) 3 > 

|<a:® fat® 

Ni-i Ng-i 

■ £ exp(-i£wrk/Ni> £ am exp <-i£trs|/Ng) , 

k*® |“® 

where k varies vertically and I varies horizontally. The inner 
summation over I is the sum over the k** 1 row for each k. The 
rule implied in this expression is that row transforms are 
taken first and then the column transforms. Let 

Ng-1 

tfks “ £ am exp (-iSirsf/Ng) 

»«0 

Then the summation over k becomes 

Ni-i 

£ #ks exp<-i£irrk/Ni) 
k«0 

This is a sum over the s*h column for each s. Th» rule implied 
by this expression is that the column transform is next taken. 
Therefore the complete rule for the two-dimensional DFT is to 


take all the row transforms and then all the c< lurnn transforms 


of this result. Bookkeeping may be a problem! Furthermore, 
the column transforms may alternat i vely be taken first and then 
the row transforms of that result. 

MATRIX APPROACH 

Following the approach of Brigham (1974), the Fourier 
transform CA] of original data Ca] may be written as 

/ A0 \ / W® U® U® U® 1 / a® \ 

< pj | = < u® U* W2 W3 |( ai I 

< As ) I W® W£ W4 W6 | I as ! 

\ As / \ W® W3 W6 / \ as / , 

with U * exp < iSir/N) . This matrix equation describes the 
transformat ion of one column vector (ai) into anr' her column 
vector (Ai>. The CU1 matrix performs the column transform by 
pre-mult i plying the column vector (ai). If the CUD matrix 
premultiplies a matrix Cal, it will transform each column of 
that matrix. The corresponding transform operator matrix which 
would give a row transform is CUT], which post-mult ipl ies a row 
vector (ai)T = (a©, ai, as, a 3 > . If CUT] post-mult i pi ies the 
matrix Cal, it transforms every row of Cal. Therefore for a 
two-dimensional matrix transformat ion of the matrix Cal, 

CUD Cal CWTl - CA1 

This shows that either row transforms or column transforms 
could be calculated first, because matrix multiplication is 
associative: 

(CU3 Cal) CUT] » CWD (Cal CUT]) . [ft] 

For a square matrix, the dimensions of CUT] are the same 
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as those for CW1 and CW^l «* CW1. If the matrix Cal is not 
square, say Ni rows and N£ columns, then the pre-mult iplying 
matrix CW1 will be Ni by Ni and the post-mult iplying matrix 
CW^l will be N£ by N£. 

The following diagram illustrates the two equivalent 
paths which may be used to obtain the two-dimensional Fourier 
transform. 

function Cal row transform -* Cal CW^i 

(post-mult i ply CW^i) 

♦ * 

column transform column transform 

(pre-mult i ply CW 1 ) (pre-mult i ply CWl) 

+ + 

CWl Ca] -* row transform -> S-D transform CWlCalCW^l 
(post-mult iply CW T 1 > 

TWO-DIMENSIONAL REAL DATA 

For real data, the transform of two rows at a time can be 
obtained using one complex transform (Brigham, 1974 ). Whether 
or not this is done, we begin by assuming the data are a matrix 
of size Ni by N£ arranged in Form Ca: 

r 30 , 0 ... a0, |\J£— 1 1 

I I 

| ■ • I 

I ■ a I 

I • ■ I 

I I 

L aNl-1,0 ... &N1-1 « NS-1 J 

Let tfij represent the matrix elements for the intermediate step 
after taking row transforms of the original matrix Cal. Then 
the one-dimensional row transforms of the real rows of Cal 
produce 
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# 

r # 0,0 # 0 , 1 ... # 0 , <N£/£)-l #0»N£/£ > #0» <N£/£)-l ••• 

I < 

I > . 

Is { * 

I > . 

| < 

I > * 

L #Ni-l»0 #N1-1, 1- » ■ #Nl~li <N£/£>~1 #Nl-l»N£/£< #N1~1* <N£/£>~1 ••• 


The redundant part is the right side with the complex 

conjugates. It is written in terms of the non-redundant part. 

The column just to the left of the redundant part, column Ng/£, 

is called the boundary element column. Using Mi = Ni/£ and Ms 

* Ns/£, the # matrix rnay be relabeled as 

* * 


r «0, 0 . . 

I • 

I ■ 

9 * 

I «Ml-lt0.-- 

I #±M1»0 .. 

I * 

I #— Ml + li 0* • 


#0. M£- 1 

#Ml-liM£-l 
#±MliM£-l 
#-Ml+l» M£— 1 


#0» ±M£ < 

> 
< 
> 
< 
> 

#M1-1» ±M£< 

> 

#±M1» ±M£> 

< 

#-Ml+l» ±M£> 


#0»M£-t ••• *0» 1 1 


« * 

#M1-1»M£-1 #Mi-l»l 

* * 

#±M1*M£-1 ••• #±M1*1 

* * 


#-Ml+liM£-l #-Ml + l»l 


I < I 

Is ) • I 

| ■ ( • I 

I s ) si 

I < * * I 

L «-l,0 ... #-iiM£-l «-l»±M£ > #-liM£-l #-l» 1 J 

Elements of the first column and the boundary element column 

are real because the original data were real. > 

Next the column transforms are taken to obtain 


# 

#0, 1 1 


* I 

#Nl-li 1 J 


Form 4ft 


# * 

r o®,0 ft®, 1 ... ft®, MS-1 ft®»±M2 < ft®, MS-1 ••• ft0* 1 

I > * # 

! < ft-1, MS-1 ft-1* 1 

i • ) • 

i . i < ii 

i > 

i < # # 

IftMl-l,® ftMl-l, 1 ••• ftMl-liMS-1 ftMl-l*±MS> ft-Ml+l,ME-l ••• ft-Ml+1* 

I < # * 

|ft±Ml,® ft±Ml,l ft±Mli MS-1 ft±Ml»±MS> ft±Ml*M2-l ft±Ml, 1 

I < * * 

I > ftMl-liMS-l fiMl-l» 1 

I < 

I . IV > III 

I < 

i > # * 

L ft-1,® ft— 1, 1 ... ft-1, MS-1 ft-1 , ±MS < ftl, MS-1 ••• ftl* 1 


The first column and the boundary element column are both 
Hermit ian. The redundant part is written in terms of the 
complex conjugates of the non-redundant part. (Note that the 
low-frequency edges of the quadrants are not at the center but 
at the far corners. ) The first subscript of the redundant part 


1 


of a row is usually not the same as that of the non-redundant 
part. The explanation is given as follows. In one dimension 
the continuous Fourier transform of f*(x) is F*(-s), and in two 
dimensions the discrete Fourier transform of #*|k is ft*-|,k. 

Except for the first row, quadrant II (III) of the Cft3 
matrix above comes from conjugating and rearranging part of 
quadrant IV (I). This symmetry can be seen readily in 
photographs of two-dimensional optical transforms, for example, 
%ri Andrews (197®), pages 34, 35, 42, 107-109, 165-167. In 
these photographs, the logarithms of the Fourier transforms ars 
shewn in order to decrease the dynamic range of the transform 
and allow the small values to be seen With the limited dynamic 
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range of an image display. 

Recall that either the row transform'- or the column 

transforms can be taken first for the two~dime‘'sional Fourier 
transform. If the column transforms are obtained first, two 
real rows (in contrast to columns) are obtained. When nex* - the 
row transforms are taken, two Hermit ian rows (in contrast to 
columns) are obtained. The final transform matrix CA3 is the 
same as that obtained by taking row transforms first and then 
column transforms, and it may therefore be written as 




Hermit ian 

Hermit ian 




co 1 urnn 

co 1 urnn 




♦ 

♦ 


Hermit ian 

4 

T (real) 

1 ( re* 1 ) > 

1 

row 


1 

< 




i 

1 ) 









1 

1 ) 


Hermit ian 

4 

1 (real) 

(real ) < 


row 


1 

1 > 




1 

< 




1 

1 > 




L 

< 

J 


redundant 

The real elements are indicated. Although the right "half" is 
shown to be redundant, the bottom "half" may be taken as the 
redundant part, as would be natural after taking column 
transforms. 

Form 4A can be rearranged into Form £A, either by 
function domain phase multiplication as previously defined, or 
equivalently by transform domain rearrangement . Transform 
domain rearrangement gives 


Form 2A ; 




# 

A±M1»MS-1 


A-l, ±M£ 
A0, ±M£ 


# 

A 1 , MS- 1 • • ■ 
# 

A0» MS-1 


L flMl-lt IMS 


# 

A-Ml + 1 » MS— 1 ••• 


# 

A±M1i 1 ^±M1*0 < ^±M1» 1 A±M1»M£-1 1 

> I 

< . I 

> . I 

< . i 


# 


> 


Al, 1 

A-l,0 

( A— 1,1 • ■ • 

A— 1 , MS— 

* 


> 


A0, 1 

A0, 0 

< A0, l ... 

A0, MS-1 


> i 

< . i 

> . i 

< . i 

# > i 

A-Ml + lil Ami- 1,0 < OMl-lil ••• AM1-1»M£-1 J 

redundant 


The dashed line is the redundant line from Form 4A. Form £A 
should be rewritten in terms of known values, so that the 
redundaxit part is given in terms of the non-red undant part. 

Form £A: 


r ft±Mlf±MS 


A-l, ±MS 
O0» ±M2 


L«M 1 - 1 » ±M£ 


A±M1»-M3+1 ••• A±M1»-1 A±M1»0 > A±M1i-1 


A±M1i -M2+1 


III 


> 

< 

> 

< 


IV 


A-l, -1 
A0, -l 


A-l, 0 
A0,0 


> * 

{ Al , — 1 ••• 
> * 

( R0i — i ... 
> 

,< 

> I 
< 

> * 

AM1-1i-M£+1 Am 1 - 1 * - 1 A|y|i-i f 0 < A-M1 + 1»-1 


A— 1 , —MS+1 ■ ■ 
A0, —MS+1 ■ • ■ 

II 


Al,-M£+1 

# 

A0, -MS+1 


• A-mi+i, -MS+ 1 J 
redundant 

If positive quadrants are desired rather than negative for the 
non-redundant part, a phase shift should be performed only on y 
for the data (or rearrangement on v in Form 4A). 
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The resulting form is: 
T A±mi» ® A±r*1l! 1 


IV 


A*M1»MS-1 A±M1,±M£ > A±M1,M£-1 

< 

> 

< 


* 

A±M1» i 


A-i» 0 
A0» 0 


A-l» 1 
A0, 1 ■ 


A-l * MS-1 
A0» MS-1 


A-l, ±M£ 

A0» ±M£ 


I 


> 

< * 

> Alt MS- 1 
< * 

> A®, MS-i 

< 

> 

< 

> 

< * 


# 

Al, l 
* 

A0, 1 


L A(y|l“l»0 AMI -1 » 1 ■■■ AMI-1 » M£— 1 AM1~1,±M£ > A-M1 + 1? M£-l A-Ml+1,1 J 

redundant 

The required phase multiplication in the function domain is 
with (-1)1, not <-l)k+I. Only the distance from the x axis 
<the y measurement) is used. This puts the high frequencies 
together and the low frequencies together in the non-red undant 
pa»?t of the transform and keeps only positive frequencies. The 
boundary-element column contains the highest frequencies. The 
phase multiplication can be done either before or after taking 
the row transforms to get the ot matrix. 




redundant 


SUMMARY PICTURES 


For complex two-dimensional data, there are two possible 
paths to the correct transform representat ion if the origin is 
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in the center. In the following diagram* the actual origin of 
the two-dimensional data is marked with CD, the origin assumad 
by tha FFT algorithm is marked with <), subscript 1 on quadrant 
numbers indicates the actual quadrants, and subscript £ 
indicates the quadrant numbers as assumad by the FFT algorithm. 



Form 4a 


Thara ara also two possibla paths for obtaining tha transform 
origin in tha centari 



It is possible that parts of both of these pictures will be 
applied to the data. 










* 


The notation of Dppenheim and Schaffer will b® used. 
Consider two array*: xi(m,n) of dimension Mi x Nj., and xg(m, n> 

of dimension Mg « Ng. The linear convolution may be defined as 

00 00 

X3(m,n) » £ E xi<q,r) xg ( (m-q> , (n-r > ) . 

q »-00 r*-00 

The limits on the summations depend on whether the convolution 
is expandirifl or non-expanding. 

It is possible to define for computational purposes 
either expanding or non-expanding convolutions. As an example 
in one dimension, consider two functions: f (of length nf) &*"<d 
g (of length ng). The convolution h « f * g will have length 

n h * nf + ng - 1, which is the normal expanding convolution 
result. For the non-expanding result, let f be the input, g 
the filter, and the output required to have the same length as 
the input, i.e. , nh “ nf. The expansion of f caused by 

convolution with g is dropped. To save computer time it would 
not be calculated at all. A non-expanding result is important 
in (1) iterative calculations, where the output would expand 
with each iteration, (2) two dimensions, when computer storage 
for only a certain size image is available, and (3) when one is 
interested in the filtered result only over the domain of the 
input. FFT programs might not conveniently accommodate the 
expanded arrays, especially in two dimensions. The limits of 
the summations in the linear convolution do not depend on the 
output variables for the expanding convolution; for non- 
expanding, they do. Note that the size of xs(m, n) defined 
above will be (Mi+Mg-1) x (Ni+Ng-1) for an expanding 
calculat ion. 
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PERIODIC CONVOLUTION 


In on« dimension, cyclic convolution was pictured with 
the two functions being convolved written on the circumference 
of cylinders which were rotated with respect to each other. 
The model for two dimensions could be one sphere inside 

another, with the spheres rotated with respect to each other. 

The periodic convolution is defined to be 
~ M— 1 N-l ~ ~ 

M3<m,n) * S £ xi<q,r) xg< <m-q> , (n-r> > 

q a 121 r»)Zi 

M and N must be specified along with the sum to completely 
define the periodic convolution, since the results of the 
periodic convolution will depend on the choice of M and N. For 
this output, m ■ 0 , 1,...,M-1 and n ■ (9,1,..., N-l. However, 

X3<m,n) need not be calculated in this order, i.e. , from 0 to 
M-l and 0 to N-l. J,t may be calculated in any order. 

If xi <m, n) has dimension Mi x Ni and X£<m, n) has 
dimension M 2 x N 2 , then the periodic convolution can be 
calculated for M £ Ml and M l M 2 and N l Ni and N It N2* The 
functions xi and X2 must be packed with zeros to make them M x 
N periodic. Therefore xi and X 2 must be filled to the same 
size to define the periodic convolution. Then the result of 
the periodic convolution will correspond to that obtained from 
the DFT Convolution Theorem, using the DFT' s of xi and X2, each 
of size M x N. xi and X2 must be packed with zeros to this 
size. Then 

3{X3<m,n>> ■. Xi<k, I) X2<k, »> 

For the periodic result to agree with the linear, M i Mi + M 2 - 
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